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Abstract. We offer a piece of evidence that the problems of finding the number of 
mutually unbiased bases (MUB) and mutually orthogonal Latin squares (MOLS) might 
not be equivalent. We study a particular procedure which has been shown to relate 
the two problems and generates complete sets of MUBs in power-of-prime dimensions 
and three MUBs in dimension six. For these cases, every square from an augmented 
set of MOLS has a corresponding MUB. We show that this no longer holds for certain 
composite dimensions. 



1. Introduction 

Mutually unbiased bases (MUB) encapsulate the concept of complementarity in 
quantum formalism. Quantum observables associated with these bases are maximally 
complementary in the sense that given a system in an eigenstate of one observable, 
measurement outcomes of the other observables are completely random. Although the 
complementarity is a distinguishing feature of quantum mechanics, we still do not know 
what the total number of mutually maximally complementary observables for a general 
(i-level system is. 

It is known that for d being a power of a prime, there are d + 1 MUBs [H [2] and 
this number sets the upper bound for arbitrary d. For all other dimensions, it is a 
puzzle whether this bound is saturated. Solving this problem gives insight not only 
into physics, but also to mathematics, as the problem is linked with other unsolved 
mathematical problems [31 H]. It was also noticed that it is similar in spirit to some 
combinatorial problems [5l El El [8] and problems in finite geometry [9l [10] • Here, we 
shall briefly review and study one such connection [8]. 

The problem which is assumed to be connected to finding the number of MUBs is 
that of finding the number of mutually orthogonal Latin squares (MOLS). The latter 
has a long history originating in the works by Euler |llj and more is known about it 
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than about the number of MUBs. For example, it is known that there are no more than 
three squares in the augmented set of MOLS of order six [12j, but the question whether 
or not there are more than three MUBs in dimension six is still open. The connection 
of Ref. [8] allows to link every square in an augmented set of MOLS with a MUB for 
power-of-prime dimensions and dimension six. Here we show that this connection fails 
in composite dimensions for which MacNeish's bound is not tight [13]. We study in 
detail the case of ci = 10, being the smallest d with this property: while there are at 
least four squares in an augmented set of MOLS, one cannot find more than three MUBs 
using the link of Ref. [8] . 

2. The connection 

A Latin square of order d is an array of numbers {0, ...,d — 1} where every row and 
every column contains each number exactly once. Two Latin squares, A = [Aij] and 
B = [Bij], are orthogonal if all ordered pairs {Aij,Bij) are distinct. There are at most 
d — 1 MOLS, and such a set of MOLS is called complete. Complete sets of MOLS are 
known to exist for d being a power of a prime. It is also known that there are no two 
MOLS in dimension six [12]. The existence of L MOLS is equivalent to the existence 
of a combinatorial design called a net with L + 2 rows [13]. The net design has a form 
of a table in which every row contains d"^ distinct numbers. They are grouped into d 
cells of d numbers each, in such a way that the numbers of any cell in a given row are 
distributed among all cells of any other row. The additional two rows of the net design 
correspond to orthogonal but not Latin squares, with the entries A^ = j and Aij = i. 
The set of all L + 2 squares is referred to as the augmented set of MOLS. An algorithm 
to construct the design from a set of MOLS is given, e.g., in Ref. p]. 

The MUBs are constructed using the entries of any cell of the design. We write the 
entries in modulo d decomposition such that each of them is now represented by two 
integers: m and n, having values from the set {0, 1, . . . ,d—l}. These integers are taken 
as exponents of Weyl-Schwinger operators, X^Z^, defined as: 

Zrfl/t) = ?7^|k), Xdln) = \{k + I) mod d) , 

with rjd = exp {i2TT / d) being a complex di\i root of unity. The Weyl-Schwinger operators 
span an unitary operator basis which is orthogonal with respect to trace scalar product. 
If one can find sets of d of them which commute (each set including identity) the joint 
eigenbases of the commuting operators form MUBs [151 [16]. K turns out that for prime d 
and for dimension six, the d operators having exponents from a single cell of the design 
commute and therefore define MUBs [8]. For power-of-prime d = p^, the operators 
having exponents from some cells do not commute. In order to obtain a complete set of 
d+1 MUBs one needs to take advantage of the fact that d can be factored. In this case, 
every integer m and n from the net design can be represented by r digits having their 
values from the set {0, 1, . . . ,p — 1}, i.e. there is a mappinjjl m i-^ (mi, m2, . . . , nir) and 

f Note that there are many such mappings. Assuming that m and n can be decomposed independently 
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n (ni,n2, . . . ,nr). We take these integers as exponents of tensor product operators 
X^^Z^^^XpZ^^ ® . . . OX^-Z;-. For suitable decompositions, related to finite fields, 
we find that the operators having exponents from a single cell of the design again 
commute and hence define MUBs. 

3. MacNeish's bound 

MacNeish gave a lower bound on the number of MOLS [13] . If two squares of order a are 
orthogonal, A B, and two squares of order b are orthogonal, C -L D, then the squares 
obtained by a direct product, of order ab, are also orthogonal, A x C -L B x D. This 
implies that the number of MOLS, C, of order d = p^^ . . .p^, Pi being prime factors of 
d, is at least C > minj(p[* — 1), where pi' — 1 is the number of MOLS of order p[\ 

A parallel result holds for MUBs [T6l [17], which we call the quantum MacNeish 
bound. If I a) and \b) are the states of two different MUBs in dimension di, and |c) and 
\d) are the states of two MUBs in dimension d2, then the tensor product bases \a) ® |c) 
and \b) ® \d) form MUBs in dimension did2. Thus, for d = p[^ . . .pj^" there are at least 
mini(p[" + 1) MUBs. 

Our motivation to study dimension ten comes from the fact that it is the simplest 
case in which MacNeish's bound is not tight. There are at least two MOLS of order ten, 
which is larger than MacNeish's bound of one. If the connection established in Ref. [8] 
holds generally, we shall correspondingly expect the quantum MacNeish bound not to 
be tight for d = 10. It is already known that the quantum bound is not tight in general, 
but the smalltest case for which it was proven is d = 26^ = 676 [7]. 

4. Ten dimensions 



The two MOLS of order ten read [T8] : 



0123456789 
1265809347 
2940573861 
3497685102 
4378160295 
5836297014 
6519738420 
7052914638 
8704321956 
9681042573 



0249187536 
1734592608 
2387641950 
3952470861 
4561928073 
5620819347 
6178304295 
7493056182 
8015763429 
9806235714 



gives (d!)^ maps. 
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where we present pairs of numbers mn. Writing the pairs as exponents of operators 
X^Z^Q, the first row gives ten commuting operators Z^q and therefore defines the 
eigenbasis of Ziq. Similarly, the second row gives the eigenbasis of Xio, and the third row 
provides the eigenbasis of XiqZiq. However, the operators corresponding to the fourth 
row do not commute, e.g. [X^qZ^q, X^qZ^q] 7^ 0, and in this way we do not improve 
upon the quantum MacNeish bound, which is three for d = 10. 

Similar to the case of d being a power of a prime, one may ask if there is a 
decomposition of m and n into pairs mi m2 and ni 71,2 (with mi,ni G {0,1} and 
m2, n2 G {0, 1, 2, 3, 4}), respectively, such that the operators X^^Z^^ (g) X^'^Z^^, having 
their exponents from the corresponding entries of the rows of ([T]), commute. However, 
contrary to the case of power-of-prime d, for d = did2 with coprime factors, the problem 
of finding commuting operators X]^^ Z^^ ^^d^^ equivalent to the problem of finding 

commuting operators X^Z^. This is a consequence of the lemma in the appendix, which 
states that the tensor product operators and the operators in d dimensions are related 
by a permutation. 

We checked by grouping the commuting operators X|^Z"q, that their eigenbases 
lead to at most three MUBs. Similarly, we verified for = 35 that there are at most six 
MUBs formed by the eigenbases of X^Zgg [19]. These are independent proofs of special 
cases of the result by Aschbacher, Childs and Wocjan [20]. In particular, they showed 
that the eigenbases of Weyl-Schwinger operators do not lead to more MUBs than given 
by the quantum MacNeish bound. Therefore, for all d in which MacNeish's bound on 
the number of MOLS is not tight, the connection of Ref. p] fails to relate all the squares 
with MUBs. 

5. Conclusions 

We presented some evidence that the physical problem of the number of MUBs might 
not be equivalent to the mathematical problem of the number of MOLS. The proof of 
this statement is still an open question. 
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Appendix 

Lemma: If d = did2 with gcd((ii, ^2) = 1; there exists a permutation matrix T such that 
fXdf-^ = Xd, ® Xd, and tzf+'^'f-' = Z^, ® Z^, . 

Proof : Define the permutatioii§| matrix T by: 

f\j) = |(j mod di)d2+j mod rfa) = \j mod di)\j mod rfa) = \ji)\j2)- 

Hence we have 

d-1 

fXaf-' = J2f\{j + 1) mod d){j\f-' 

j=0 

di-l d2-l 

= \Ul + 1) mod di){ji\ ® I (-^2 + 1) "^Od d2){32\ 

jl=0 j2=0 

= (g) Xd2 , 



and 



-1 

^^d,+d2f-i ^ J2Vd'^'''^'T\j){j\f-^ 

j=0 

d-1 



j=0 

di-1 d2-l 
$^<|jl)(jl|®$^^£|j2)(j2| 

ii=o j2=o 

Zdi ® Zd2, 



where we used rj^, = rj-jj^^"'^ = 77 j. for i = 1, 2. 

References 

[1] Ivanovic I D 1981 J. Phys. A 14 3241 

[2] Wootters W K and Fields B D 1989 Ann. Phys. 191 363 

[3] Boykin P O, Sitharam M, Tiep P H and Wocjan P 2005|arXiv: quan t-ph/0506089| 

[4] Bengtsson I, Bruzda W, Ericsson A, Larsson J-A, Tadej W and Zyczkowski K 2007 J. Math. Phys. 
48 052106 

[5] Zauner G 1999 Dissertation (Universitat Wien) 

[6] Bengtsson I 2004 |arXlv:quanT^'h/M06174| 

[7] Wocjan P and Beth T 2005 Quant. Inf. Comp. 5 93 

[8] Paterek T, Dakic B and Brukner C 2009 Phys. Rev. A 79 012109 

[9] Saniga M, Planat M and Rosu H 2004 J. Opt. B 6 L19 

[10] Wootters W K 2006 Found. Phys. 36 112 

[11] Euler L 1849 Comm. Arithm. 2 593 

[12] Tarry G 1900 Comptes Rendu de I'Assoc. Francaise pour V Advancement de Sci. Nat. 1 122 



§ It follows from the Chinese remainder theorem that T is a permutation, and this is why d has to 
have coprime factors. 



On the connection between mutually unbiased bases and orthogonal Latin squares 6 
[13] MacNeish H F 1921 Ann. Math 23 221 

[14] Colbourn C J and Dinitz H J 1996 The CRC Handbook of Combinatorial Designs (CRC Press, 
Boca Raton, FL) 

[15] Bandyopadhyay S, Boykin P O, Roychowdhury V and Vatan F 2002 Algorithmica 34 512 
[16] Grassl M 2004 arXiv:quant-ph/0406175 

[17] Klappenecker A and Rotteler M 2004 Lect. Notes Comput. Sci. 2948 137 
[18] Parker E T 1959 Proc. Nat. Acad. Sci. 45 859 
[19] Wojtas M 1996 J. Comb. Des. 4 153 

[20] Aschbacher M, Childs A M and Wocjan P 2007 J. Alg. Comb. 25 111 



